Accelerator modes of square well system 
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We study accelerator modes of a particle, confined in an one-dimensional infinite square well 
potential, subjected to a time-periodic pulsed field. Dynamics of such a particle can be described 
by one generalization of the kicked rotor. In comparison with the kicked rotor, this generalization 
is shown to have a much larger parametric space for existence of the modes. Using this freedom 
we provide evidence that accelerator mode assisted anomalous transport is greatly enhanced when 
low order resonances are exposed at the border of chaos. We also present signature of the enhanced 
transport in the quantum domain. 



I. INTRODUCTION 

In recent years there is a surge of interest in studying 
mixed systems, with both regular and chaotic regions co- 
existing in phase space. Mixed systems are important as 
a majority of dynamical systems found in nature fall in 
this category. One of the outstanding problems is to re- 
solve how chaotic orbits are influenced by regular regions 
in the long time limit. In this direction of study, two 
dimensional area-preserving maps are promising candi- 
dates as they exhibit many generic features. A text book 
paradigm of this kind is the standard map of kicked ro- 
tor 0]. For certain values of system parameter, this map 
possesses special kind of orbits called Accelerator Modes 
(AM) (2] . If an initial state of the rotor corresponds to 
the AM, then it is uniformly accelerated after each kick. 
In other words, momentum of the rotor increases linearly 
with time. This is to be contrasted with generic chaotic 
orbits that diffuse in momentum, i.e., momentum vari- 
ance of an ensemble of chaotic orbits increases linearly in 
time leading to normal diffusion. 

AM are seen in phase space as small regular islands 
embedded in a chaotic sea. They are separated from 
chaotic region by cantori (partial-barriers). If a chaotic 
orbit enters through the cantori, it sticks to the cantori 
for very long time and thus is accelerated along with the 
modes. Considering n as time, such a longexcursion of a 
chaotic orbit is an example of Levy flight Levy flight, 
some times called as super diffusion, is characterized by 
(Pn) ^ "-^ with 7 > 1. AM assisted anomalous transport 
of this kind has been already studied with different moti- 
vations 4] . It should be noted that the above mechanism 
of anomalous effect is numerically ascertained. 

Further studies have shown that these modes influence 
corresponding quantum system as well 5, 6J. In partic- 
ular, the existence of regular islands in chaotic sea sig- 
nificantly alter the localization of quasienergy states in 
angular momentum, leading to enhancement in quantum 
transport. Developments on confining atoms in magneto- 
optical traps have opened experimental feasibility for 



atom-optics realizations of the kicked rotor. One such 
realization, in which an ensemble of cold cesium atoms 
are exposed to a periodically pulsed standing wave light, 
has shown that for certain values of the kick amplitude 
the momentum distribution of the atomic sample is non- 
exponentially localized This then is shown as an ef- 
fect of AM in the quantum system, in contrast to the 
generic exponential localization. 

Semiconductor based technological developments, on 
the other hand, have shown that it is possible to fabri- 
cate potential wells on atomic scales. Motion of electrons 
in such quantum wells in presence of an external electro- 
magnetic field has been studied for experimental signa- 
tures of quantum chaos Q. A simple model of this kind 
is a particle confined in a one-dimensional infinite square 
well exposed to a time periodic pulse field. This system 
is known to be one generalization of the kicked rotor . 
The virtue of this generalization arises from two length 
scales of the system, namely, width of the well and field 
wavelength. If the well width is an integer multiple of 
the field wavelength, kick to kick dynamics of the par- 
ticle is equivalent to that of the rotor. Otherwise, the 
dynamics is described by a discontinuous map resulting 
in a new scenario for transition to chaos even for weak 
field strengths 0. 

In this paper we explore this model in the context of 
AM. We show that, our generalization leads to a larger 
parametric space, in comparison to the standard map, 
for existence of the AM. We demonstrate that when the 
border of AM (which we call a "beach" ) exposes a chain 
of islands of a lower order resonance, "stickiness" of the 
beach increases significantly leading to anomalous trans- 
port. We also study the signature of such classical phe- 
nomena in the corresponding quantum system. 



II. CLASSICAL SYSTEM 

A. The Map 

Let us consider a particle confined in an one- 
dimensional infinite square well potential Vq{x) of width 
2a and the Hamiltonian is 
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The particle is subjected to a time periodic pulsed field 
of period T . The perturbed system is governed by the 
Hamiltonian 



H 



Ho + ecos (^^^^6 



(2) 



where e and A are strength and wavelength of the field 
respectively; a train of delta functions accomplishes the 
time periodic pulse. Kick to kick dynamics of the particle 
can be described by the map 



Xn+i = (-lf"(X„ + /^„) + (-l)""^^Sgn(P„)B 



-Bn + lt 



Pn+1 = (-lf"/^„ + ^sin(2^i?X„+i), (3) 

where B„ = [Sgn(P„)(X„ + P„) + 1/2]; [..] and Sgn(..) 
stand for integer part and sign of the argument respec- 
tively. Bn is the number of bounces of the particle 
between the walls during the interval between nth and 
(n-l-l)th kick. This dimensionlcss map is related to phys- 
ical variables by the scaling: 

where K is effective field strength and R is ratio of two 
length scales of the system. Henceforth we refer to the 
classical mapping (O as the well map. 

The well map can be studied by invoking a Generalized 
Standard Map (GSM) 

K 



K 



Jn+i = Jn + sin(27rP6'„) , 
6*71+1 = dn + Jn+1 (mod 1) . 



(5) 



since time reversal of GSM and ((SJ are quantitatively 
related. GSM is defined on a cylinder (— oo, oo) x 
[—1/2, 1/2) and the standard map is a special case with 
P = 1. A detailed study of GSM can be found in our 
earlier work The GSM is highly chaotic for strong 
field strength (large K). In addition, it exhibits chaotic 
motion even for weak field strength (small K) depending 
on the parameter R. In order to study the AM in the 
well map, we invoke the GSM which is amenable to a 
detailed analysis. 



B. Stability Analysis 

Periodicity of the GSM in J and 9 (with unit period) 
implies existence of AM. These are located at (J', 9'): 



K 

f = — sm{2nR9') = I , 



(6) 



where m and / are integers. They are also termed as step- 
1/1 AM as the acceleration is 1/1 at each iteration. Fixed 



points belong to the family of AM with / = 0. For stable 
AM, required stability condition is 



KR 

and using © this becomes 



< cos(27rP6l') < , 



(7) 



(8) 



This inequality can also be rewritten as 

< P < Pi , (9) 

where 

2 ^ -1/2 




For stable AM 



\9'\ 



2P 



2j 



1 =F — sin 

TT 



27r|/| 



(10) 



where j is an integer. Thus the AM are characterized 
by two integers / and j and we call them as {l,j)^ type 
modes. Since the cylindrical phase space of GSM has 
the constraint \9'\ < 1/2 (equivalent to restriction of the 
particle dynamics to between walls of the well) , zero can 
not be lower limit in the above inequality. 

The parameter P in GSM is recognized as frequency 
of the sinusoidal force term. As P increases from zero 
\9'\ decreases such that type modes first appear in 

the phase space followed by the (/, j)_|_ type. In a similar 
way, as P decreases \9'\ increases such that {l,j)^ type 
is the first to disappear from the phase space while the 
type mode is the last to disappear. Hence, lowest 
limit for the inequality © is set by the disappearance 
of (/, j)„ type modes. From Eqn. Hll)|l . the constraint 
\9'\ < 1/2 for {l,j)_ becomes Pq < P where 



1 1 /27r|/| 
Po = 2i-f l--sin-M — ^ 

TT \ 

Thus the inequality © is replaced by 
Pq P < Pi . 



(11) 



We emphasize that both the inequalities © and l(TT|) 
must be simultaneously satisfied for existence of AM in 
the phase space. This is in contrast to the standard map 
for which the inequality © with P = 1 itself is suffi- 
cient. Thus the two control parameters, K and P, pro- 
vide a larger parametric space, in comparison to that of 
the standard map, for the existence of AM. At PT = 2tt\1\, 
the lower bound of the inequality ©, Po = 1/2 and this 
is the minimum possible value of Pq. In other words, for 
P < 1/2 AM do not exist. In fact, GSM is hyperbohc for 
P < 1/2 and there are only unstable orbits in the phase 
space. 
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FIG. 1; Stable regions embedded in chaotic sea are the accelerator modes of the well map (Kjl-K = 1.05) whose P' = and X' 
is positive. For K = 0.7, 0.8 and 1, we observe (1, 0)_ type mode. For K = 1.4, one more mode just appears in the phase space 
which is (1,0)^ type. For further increase of i?, both the types are fully visible. Since X' oc 1/_R, as J? increases the accelerator 
islands move towards interior of the square well. Note the three different scales on P-axis. 



Results 



Since the well map and the GSM are quantitatively 
related, location of AM are same for both the maps i.e., 
{X' , P') = (6*', J') and all the above arguments hold for 
well map also. However, AM of the well map differ only 
in the following way due to the reflective boundary con- 
dition. When particle is on the AM, nth kick causes the 
particle to undergo n bounces between the walls. If n 
is odd momentum changes its sign, in turn sign of the 
particle position is also flipped. 

Phase space of a mixed system generally has many in- 
tricate structures. In particular, boundaries of regular re- 



gion embedded in chaotic sea have structures at all scales 
and all of them are hard to discern. Stable AM islands 
are believed to be separated from the remaining chaotic 
region by cantori (partial-barriers). Phase space in the 
vicinity of the cantori is sticky as it retains long time 
correlations. When chaotic orbits explore all parts of the 
phase space, there are occasions that they pass through 
the cantori. When such an event occurs, chaotic orbits 
stick to the vicinity of cantori for long time. Thus they 
are dragged along the modes balhstically 0. This is a 
mechanism which is known to enhance momentum trans- 
port. Below we demonstrate that the beaches around the 
AM island when occupied by a large resonance of small 
order are highly sticky. This in turn implies that the 
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cantori surrounding these must have large gaps to allow 
larger penetration of itinerant chaotic orbits. 
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FIG. 2: Exponent 7 is plotted with respect to the parameter R 
for K/2ti — 1.05. To obtain the exponent, an initial ensemble 
of 10000 phase space points is iterated by the well map for 
5000 time steps. For (i) the ensemble consists of points with 
Po = and Xq distributed uniformly between —1/2 and 1/2. 
The range of R given by the Eqn. lllll with Ro = 0.598 and 
Ri — 1.988 is seen as 7 ~ 2 due to ballistic evolution of the 
accelerating points contained in the ensemble. For (ii) the 
ensemble is chosen from chaotic region of the phase space. 
Maximum exponent observed for R = 0.8, marked by circle, 
is due to the beaches around the AM island. 

Before we dwell more on AM assisted super diffusion, 
it is instructive to look at some of the modes of the well 
map. Fig. ^ shows islands for different R values. In 
the given range of phase space we find one (1,0)_ type 
mode for R = 0.7,0.8 and 1. As R increases the islands 
move towards interior of the well. For R= 1.4, one more 
mode just appears at boundary of the well. This is a 
(1,0)_,_ type mode. For further large values of R both 
the types fully appear and they move into the interior of 
the well. Evidently, parameter R also changes the size of 
the islands. 

In order to quantify the role of R in AM assisted 
anomalous transport, we consider ((P„ — Pq)^) ~ n"' . 
Here the angular bracket stands for an ensemble aver- 
age and the exponent can be evaluated by iterating an 
ensemble of phase space points for long time under the 
well map. We perform numerical calculations for two dif- 
ferent initial ensembles: (i) phase space points uniformly 
distributed along the line Pq — which correspond to 
both accelerator and chaotic orbits; (ii) phase space points 
from a chaotic region which include only chaotic orbits. 

The ensemble (i) contains more chaotic orbits than the 



regular (accelerator) orbits. Chaotic orbits exhibit nor- 
mal diffusion in momentum i.e., like the random walk, 
until they are dragged along the mode. On the other 
hand, the accelerator orbits exhibit ballistic behaviour 
(7 = 2). Hence, in the long time limit, evolution of ac- 
celerator orbits dominate the ensemble average. Fig.|2Ii), 
shows the characteristic exponent 7 for different R val- 
ues. As we see, for the range of R in Eqn. (|ll|l 7 ~ 2, 
confirming the existence of AM islands. Here the lesser 7 
for i? = 1.5 is attributed to small size of AM (see Fig.^. 
On contrary, for parameters outside the range of Eqn. 
(|ll|l . diffusion is normal as there are no AM in the phase 
space. 

Fig. Olii) corresponds to the ensemble which contains 
only chaotic orbits. Within the parametric range where 
AM exist, the exponent 7 shows large fluctuations. There 
are also occasions wherein presence of AM does not en- 
hance the diffusion significantly. It is natural to attribute 
the observed fluctuations to size of the islands and stick- 
iness of their neighbourhood. That is, big AM islands 
with very sticky neighbourhood can significantly enhance 
the transport. 

For R — 0.8, the exponent 7 is found to be maximum 
(~ 1.6). As seen from Fig.^ in this case AM island is sur- 
rounded by a chain of 1/5 resonance zones. These zones 
are also step-1 accelerators and they are separated from 
chaotic region by their boundary, the so called "beaches" . 
Dynamics on the beach is chaotic, but very sticky with 
longer classical staying time. If a wandering chaotic or- 
bit happens to approach the beach region, it stays there 
for long time and thus gets accelerated along the modes 
resulting in very large enhancement in momentum trans- 
port. This maximum transport is a generic behaviour 
i.e., independent of choice of the chaotic ensemble. From 
Fig.^we note that, although size of the AM for R = 0.7 is 
larger than that for R — 0.8, exponent 7 is maximum for 
the latter. This shows that the emergence of lower order 
resonance zones are chiefly responsible for the stickiness 
of the beach regions, leading to AM assisted anomalous 
transport. Thus the well map or equivalently GSM has a 
control parameter R by tuning which AM induced super 
diffusion can be maximized. 

More over, longer the period of resonance zones lesser 
the stickiness of associated beach regions. Upon magni- 
fying boundaries of AM islands for R — 0.7 and 1, we 
found long chain of resonance zones and their beach re- 
gions. As seen from Fig. [5|(ii), enhancement is limited 
for both these cases. 

Further, we explore step-2 AM which are smaller in 
size than the step-1 modes. We choose K/2tt — 2.05 and 
step-2 modes exist in the range 0.57 < R < 1.41. Nearly 
ballistic evolution in this range and the otherwise normal 
diffusion are shown in Fig.l^i). In Fig.lJfii) the exponent 
exhibits large fluctuations. In order to appreciate the 
role of beach regions in large scale transport, we focus 
our attention on two distinct cases where the exponent 
shows a minimum and a maximum. 

As we see from Fig. 0] when the enhancement is mini- 
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FIG. 3: Exponent 7 is plotted with respect to the parameter 
R for K/2n = 2.05. Rest of the calculation details are similar 
to the Fig. 121 In (i) the range of R given by the Eqn. Ulljl 
with _Ro = 0.57 and Ri — 1.41 is seen as 7 ~ 2. In (ii) 
the exponents for R = 0.8 and 0.9 are marked by circles. 
Corresponding AM islands for these two cases are shown in 
Fig.Bl 



mum (or absent) AM islands do not have the surrounding 
beach regions. On the other hand, transport is maxi- 
mum when AM island is surrounded by short chain of 
resonance zones and their beach regions. These results 
reinforce that the expose of lower order resonance zones 
around the AM are chiefly responsible for stickiness of 
the beach and thus the resulting enhanced momentum 
transport. 



III. QUANTUM SYSTEM 
A. The Map 

Analogous to the classical system, the kick to kick 
quantum dynamics can be studied using quantum map 
by solving Schrodinger equation for one period T as 



(12) 



where U = exp{—ikcos{2TTx/X)}exp{—iHoT/h} with 
k — eT/h, being the effective field strength. Eigenvalue 
equation of the unperturbed system is i?o|0n) = En\4>n)- 
The eigenstates and eigenvalues are given by 



a/2 cos(7i7rX), for n odd 
\/2sm(mrX), for n even 
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FIG. 4: Step-2 AM are shown for K/2n = 2.05. For R = 0.9, 
the AM island is accompanied by a chain of 2/7 resonance 
zones and sticky beaches. On the other hand, AM island for 
R = 0.8 is not accompanied by any beach regions. 



En 



2M 



(13) 



with n = 1, 2, 3 . . . Considering unperturbed eigenstates 
as the basis, time evolution of an arbitrary quantum state 
is l^Pit)) = EnMtMn) where A^it) = (0„|V'(O>- With 
following substitution 



^ ^ 'n^ ^ 2M 
quantum map 112|) takes the form 



(14) 



(15) 



where Umn = e~''^"'((/)m| exp{— jfccos(27ri?Ar)}|0„). 

Dimensionless quantum parameters fc,T and the clas- 
sical parameters are related through K/R = Skr. Since 
k l/Ti and t ^ h, semiclassical limit for given classical 
system can be achieved by taking the limits fc — > 00 and 
T 0. In the following numerical calculations, quantum 
map H15() is implemented with a finite, say N, number of 
unperturbed basis states. By obtaining the limit 



lim 



N 
n=l 



(16) 



we ensure that such a truncation of basis does not influ- 
ence the quantum system. 

As seen from Eqn. (|15|l . unperturbed motion of the 
particle between successive kicks just adds phase to the 
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wave function components through the parameter r (mod 
27r). If T — 2TT, unperturbed motion is absent and this is 
called quantum resonance. Under this condition, without 
loss of generality, we can write 



m)) = 



'l^(O)), 



(17) 



where t represents the number of kicks. Then in the limit 
t ^ oo, kinetic energy of the particle grows quadratically 
i.e., {E)^ ~ This phenomenon is similar to that 

of the quantum kicked rotor . An earlier study has 
shown that quantum resonance of kicked rotor occurs 
when T is rational multiples of 27r |13] . This non-generic 
pure quantum phenomenon does not have any classical 
analogue. In order to investigate signatures of AM in 
the quantum system, it is necessary to suppress the res- 
onance effect. This can be accomplished by taking t as 
an irrational multiple of 2tt. 



B. The wave packet 

One natural way to study the quantum signatures of 
the AM is using the time evolution of an initial wave 
packet localized in chaotic region of the classical phase 
space. We take a Gaussian wave packet confined in the 
square well as the initial state. In position representation 
it is given by 



(X|^(0)) = Cexp<^- 



2a2 



'<^>-^j, (18) 



where cr measures width of the wave packet, centred at 
{X) with momentum (P). Here he = 2r/7r^ is the ef- 
fective Planck constant. Normalization constant, C, is 
obtained from the condition 



1/2 



c 



-1/2 



2/(70F 



\{X\m)\' dX ^ 1 



erf(y+) - erf(?/_) 



with y± 



±1/2- (X) 



(19) 

As the quantum map is described in unperturbed basis 
states, the initial wave packet is represented in this basis. 
For evaluating An(0) we consider following integral 

.1/2 

Gn^V2 e*"''-^ (^IV'(O)) dX . (20) 

J -1/2 

With change of variable u = [X ~ {X) ) / ^/2a the integral 
becomes 



Gr, - 2aC e*^"<^> 



-(^^-^^<^^r.u) J^^^ (21) 



where Zn = nir + 
standard integral 



P)/he and u± = y±l\p2. Using the 



l^/Zexpf^— erff^^h (22) 



2 V a, 



Gn is represented in terms of complex error function, 
which can be computed using the algorithm 0|. The 
wave packet in the unperturbed basis is then given by 



A„(0) 



(G„+G_„)/2, 
(G„ - G_„)/2z, 



for n odd 
for n even 



(23) 



It should noted that because of the spatial confine- 
ment, the initial state is a truncated Gaussian and hence 
it is not a minimum uncertainty wave packet. However, 
it can be made close to the minimum uncertainty wave 
packet provided mean {X) is not close to boundary of 
the well and AX <C 1, i.e., spread in position is much 
less than the width of the square well. With this con- 
dition on the wave packet we consider the uncertainty 
relation AXAP = ?ie/2 for the following exercise. We 
take AX = 0.02 throughout. 



C. Results 

The initial wave packet is chosen with {{X),{P)) ~ 
(0, 0.5) such that it is placed in a chaotic region. We 
take Tie = 0-008 and hence AP = 0.2. For these choice of 
parameters, the initial state is squeezed in position and 
elongated in momentum. The number of basis states 
considered is as high as iV = 7500. Dimensionless kinetic 
energy of the time evolved state is then given by 

(E), = {m\p'\m) 



2 2 

n . 



(24) 



which is equivalent to the classical energy (Pf). 

Shown in Fig.jSjare typical typical quantum evolutions 
for K/2n — 1.05 with varying R values. In all the cases, 
initially the quantum energy increases and then attains 
a quasiperiodic saturation. However, the saturation is 
much higher for R = 0.8 than for R = 0.7 and 1. It 
should be noted that even though effective field strength 
k is larger for R = 0.7 than for R = 0.8, energy sat- 
uration for the later is much higher. We may compare 
these results with the corresponding phase space shown 
in Fig. n Thus we obtain a clear evidence that the ex- 
pose of lower order resonance zones and the associated 
sticky beach around the AM are mainly responsible for 
AM assisted enhanced transport in quantum system also. 
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FIG. 5: Kinetic energy of the quantum particle for K/2ti — 
1.05. For 7? = 0.7,0.8 and 1 the effective field strengths are 
k = 30.3,26.5 and 21.2 respectively. 

IV. DISCUSSION 

We studied accelerator modes for a particle inside a 
ID infinite square well in presence of a pulsed external 
field. Influence of modes in altering large scale trans- 
port of the classical system has been studied by invoking 
an equivalent dynamical model, the generalized standard 
map. The virtue of the generalization arises from the ra- 
tio of two length scales, namely, well width and field wave 



length. This generalization is found to have larger para- 
metric space for existence of AM. 

We demonstrated that when the border of AM (beach) 
exposes short periodic resonance zones, stickiness of the 
beach region is significantly increased, leading to AM as- 
sisted super diffusion in momentum. It also appears that 
longer the period of resonance lesser is the stickiness 
of the beach, results to less enhancement in transport. 
The corresponding quantum system is studied using time 
evolved wave packet. We showed that the presence of 
lower order resonance around the modes are responsible 
in enhancing the quantum transport as well. 

Some investigations on standard map show that 
for certain "magic" values of system parameter, AM 
islands are accompanied with hierarchy of self-similar 
structures. The vicinity of these structures are known 
to be regions of stickiness for the chaotic orbits 0, 0| , 
causing enhancement in transport. Further, various 
multi-island structures in the near threshold regime of 
the mapare identified as trapping zones for wandering 
orbits [13 leading to Levy-like flights. All these results 
collectively support that island-around-island structures 
are greatly responsible for anomalous transport in the 
mixed phase space. Using fractional kinectic equation 
[TgL , it will be of interest to enquire the connection 
between the period of boundary islands and stickiness of 
the associated beach region. 

Acknowledgement 

One of us, VBS, would like to acknowledge the hospi- 
tality of the Institute of Theoretical Atomic and Molecu- 
lar Physics, Harvard University, Cambridge, USA, during 
August 2002 when this paper was being completed and 
written. We gratefully acknowledge very useful discus- 
sions with Dr. A. Lakshminarayan. 



[1] A.J. Lichtenberg and M.A. Lieberman, Regular and 
Chaotic Dynamics, Springer- Verlag, New York, 1992. 

[2] B.V. Chirikov, Phys. Rep. 52, 263 (1979). 

[3] M.F. Shlesinger, G.M. Zaslasvsky and U. Frish (Eds.), 
Levy Flights and Related Topics in Physics, Springer- 
Verlag, 1995. 

[4] Y.H. Ichikawa, T. Kamimura and T. Hatori, Physica 29 

D, 247 (1987); R. Ishizaki, T. Horita, T. Kobayashi and 

H. Mori, Prog. Theor. Phys. 85 (5), 1013 (1991). 
[5] J.D. Hanson, E. Ott and T.M. Antonsen Jr., Phys. Rev. 

A 29, 819 (1984). 
[6] Bala Sundaram and G.M. Zaslavsky, Phys. Rev. E 59, 

7231 (1999); A. lomin and G.M. Zaslavsky, Chaos 10 

(1), 147 (2000). 
[7] E.G. Klappauf, W.H. Oskay, D.A. Steck and M.G. 

Raizen, Phys. Rev. Lett. 81, 4044 (1998). 
[8] PB. Wilkinson, T.M. Fromhold, L. Eaves, F.W. Sheard, 

N. Miura and T. Takamasu, Nature 380, 608 (1996). 
[9] R. Sankaranarayanan, A. Lakshminarayan and V.B. She- 

orey, Phys. Lett. A 279, 313 (2001). 



[10] M.F. Shlesinger, G.M. Zaslavsky and J. Klafter, Nature 
363, 31 (1993). 

[11] R. Sankaranarayanan, A. Lakshminarayan and V.B. She- 
orey, Phys. Rev. E 64, 046210 (2001). 

[12] G. Casati, B.V. Chirikov, J. Ford and P.M. Izrailev, Lec- 
ture Notes in Physics Vol.93, 334 (1979) 

[13] P.M. Izrailev and D.L. Shepelyanskii, Theor. Math. Phys. 
43, 553 (1980). 

[14] M. Abramowitz and LA. Stegun (Ed.), Handbook of 
Mathematical Functions, p. 303, Dover, New York, 1965. 

[15] W. Gautschi, SIAM J. Numer. Anal. 7(1), 187 (1970). 

[16] G.M. Zaslavsky, M. Edelman and B.A. Niyazov, Chaos 
7 (1), 159 (1997). 

[17] S. Benkadda, S. Kassibrakis, R. White and G. Zaslavsky, 
Phys. Rev. E 55, 4909 (1997). 

[18] R.B. White, S. Benkadda, S. Kassibrakis and G.M. Za- 
slavsky, Chaos 8 (4), 757 (1998). 

[19] G.M. Zaslavsky, Phys. Rep. 371, 461 (2002). 



